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Direct CP violation in decays such as B~ — t D°K~ is sensitive to the CKM angle 03 because these decays allow the interference 
of b-quark to c-quark with fe-quark to u-quark transitions. Indeed, 4>3 may be determined if one can infer the strong phase of the 
B and subsequent D° decays from experimental data. In this talk, I will discuss how this can be carried out using either a single 
decay mode of the D° by combining data from a number of D° decay modes as well as the use of other, analogous decays and the 
prospects of implementing such methods at various B-factories. Since the properties of the D° decays are crucial to these methods, 
it is possible that D°-D° mixing at the 1% level will contaminate the results. I will therefore discuss various methods to remove 
such confounding effects so that <j>3 may be determined even if such mixing is present. 
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1 Introduction 

The asymmetric B experiments BaBaR0 and BELLE i 
have already obtained preliminary measurements of the 
angle <f>\ of the unitarity triangle U of the Cabibbo 
Kobayashi-MaskawaQ (CKM) matrix through the "gold- 
plated" modetj ipK. Using B°B° oscillation it may also 
be .-possible to extract 02 via modes such as nn and 
irpu however to extract 03 using oscillations requires B s - 
mesons and so is inaccessible to T(45) machines. 

Although experiments to extract 03 via B s oscilla- 
tions may be performed at hadronic S-facilities it is also 
possible to measure 03 through direct CP violation in 
the B system. Thus, the complete set of unitarity an- 
gles may, in principle, be accessible at T(4S r ) machines. 
Specifying as many parameters of the unitarity trian- 
gle as possible is, of course an important check of the 
Standard Model (SM). In addition, if 03 is measured via 
direct CP violation the comparison to the measurement 
through indirect CP-violation in the B s system provides 
another non-trivial check of the SM. 

The idea behind the measurement of 03 through di- 
rect CP violation is to consider a process which allows 
interference of the quark level processes b — > ucs and 
b — > ucs. This may be accomplished if botl^pracfisses 
ultimately hadronize to a common final stateLmtffiJ. In 
particular b —> ucs can drive the decay B~ — > D°K~ 
while b — ► ucs can drive the decay B~ —> D°K~ which 
will thus interfere provided that D° and D° are detected 
through decays into a common final state. 

In this talk, I will discuss various strategies for the 
determination of 03 in such decays. The most crucial 
element is the selection of the D° decay modes which 
are to be used. In the simplest case is where CP viola- 
tion is seen in a single mode we will see that there is not 
enough information to precisely determine 03. However, 
in Section 2 I will show if the CP violation is large, signif- 



icant bounds may be placed on 03 thus modes where CP 
violation may be large are especially significant. More 
generally, as discussed in Section 3, 03 may be extracted 
if two or more modes are measured. A single three body 
mode may also give equivalent information because each 
point on the D-decay Dalitz plot may be regarded as a 
separate "mode" . These methods are subject to possible 
contamination from DD oscillation if it is on the order 
of 1%, I will discuss the impact of this possibility and 
methods to deal with it in Section 4 and in Section 5 I 
will give my conclusions. 

2 One D° Decay Mode 

Consider the case where D° and 5° can decay to a com- 
mon final state X. This may either be CP eigenstate (e.g. 
K + K~) or a state such as K + ir~ which is a Cabibbo 
allowed (CA) decay of the D° but a doubly Cabibbo 
suppressed (DCS) decay of the D°. In either case, one 
can determine 03 if one can measure □ the rates d(X) = 
Br(B~ -> K~[X]) and d(X) = Br(B+ -> K+[X\) 
(here [X] means a decay to X via D° mixed with the 
D channel) provided one also knows the branching ra- 
tios a = Br(B- -> K~D°); b = Br(B~ -> K~ D°); 
c(X) = Br(D° -> X) and c(X) = Br(D° -> X). This 
information allows us to solve (up to an eight fold ambi- 
guity) for the weak phase 03 as well as the total strong 
phase difference £. 

In practice, however, b is not easy to determined di- 
rectly because it is difficult to find a prominent tag for 
D°. For instance, a leptonic tag has a large background 
from leptonic decay of B~ while if one tries to tag it 
through decays such as D° — > K + tt~ the signal is sub- 
ject to interference effects from D° — ► K + tt~. In fact 
it is the existence of such interference effects we wish to 
exploit as a source of CP violation that make the direct 
determination of b via a hadronic tag impossible. 
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Of course if CP violation were seen (i.e. d^ d) then 
within the SM it must be the case that ^ ^ 0. As might 
therefore be expected, in the absence of b we can use the 
rest of the information to establish a lower bound on 03 . 
In particular, if we define Q = sin 2 03 then we obtain the 
following bound on Q and incidentally b: 

Q > Qrmn = (1 + z) (l - y/l-y/(l + z)^ /2 

-y/l + Z + \y\ < < y/l + z-\y\ (1) 

where 1 + z = (d(X) + d(X)) / (2ac(X)) , y = (d(X) - 
d(X))/(2ac(X)) and u = bc(X)/ac(X). 

In the case where z is negative, an upper bound llil 
can also be placed on Q: 

Q<l + z (2) 

which is similar in form to the upper bound on Q ob- 
tained from tree-penguin interference in B — ► Ktt tn. 
Note that it can apply even if CP violation is absent 
however z may only be negative if it happens that u < 4. 

Motivation for these bounds can be found in Fig. [j] 
where we plot the relation between u and 03 for the ex- 
perimental inputs z = 1.5 and y = 0, 1 and 2 where the 
larger values of y correspond to the smaller "lazy eight" 
curves. The boxes indicate the bounds established by the 
inequality eq. (|l|). Since y is proportional to CP viola- 
tion, it is clear that the most stringent inequality bounds 
obtain where CP violation is large. A useful property 
of these curves is that even though the strong phase dif- 
ference is not explicitly given, it may be read off (up to 
a four fold ambiguity) since for a given value of (03 , u) 
on the curve, the horizontal line through that point also 
intersects the curve at (£',it) where one of {£' ,tt — 
7r + £'} is the strong phase difference. 

Large CP violation is only possible when the two 
interfering amplitudes are similar in magnitude. Such 
a situation may happen if we consider a final state X 
where D° -> X is DCS (e.g. K+Tr~). Thus, while a 
is about two orders of magnitude greater than b due to 
color suppression, this if offset by c(X) being about two 
orders of magnitude greater than c(x) so bc(X) ~ ac{X). 
It is obviously advantageous to experimentally study all 
modes of this kind in order to find the one which gives 
the largest lower bound on Q. 

Of course we would like to determine 03 exactly 
rather than just establishing a bound for it. There are 
three possible approaches to obtaining this quantity if 
b cannot be experimentally measured: First of all, one 
could use a theoretical model to estimate b. Second, one 
could use an analogous decay B° — > K°[D — ► K + tt~] 
or At — > A[D° — > K + ir~] where the cross channel is is 
not color enhanced and interference are ~ 30% (we will 
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Figure 1: Each of the solid lines shows the locus of points in cj>s 
versus u of allowed solutions given z; = 1.5 for yi = (outer curve), 
1 (intermediate curve) and 2 (inner curve). The boxes indicate the 
inequalities Eqn. (Ilj). Note that u oc b. 

discuss this more below) or third of all, one can consider 
multiple D decay modes each of which may decay with 
a different strong phase. In this last case I also include 
D decays to a multi-body final state where each point in 
phase space may be considered as a separate "mode" . 

3 Two Modes and Three Body Decays of D° 

If d and d are measured for exactly two modes, then 
there are the same number of equations as unknowns 
and Q may be determined up to some discrete ambiguity. 
Graphically, two general curves such as in Fig. [I] will 
intersect at up to 4 points so in this case there may be a 
4- fold ambiguity in Q therefore a 16-fold ambiguity in fa. 
If three or more modes are considered then the curves in 
the b — 03 plane should only intersect at one point in the 
first quadrant of 03 so 03 has a 4-fold ambiguity. 

Ip_Fig. U we show the results of a sample calculation 
fromEj where the modes: K + tt~ , Kstt , K + p~ , K + a^ , 
Ksp~ , K* + n~ were considered. For the parameters of 
that calculation, the inner edge of the shaded region indi- 
cates the 68% CL and the outer region indicates the 95% 
CL given N — (number of B^) (acceptance) = 10 s . In 
this example it was found that with TV = 10 8 statistical 
errors in 03 were ~ 5° — 10° for a variety of initial values 
of 03 and strong phases. 

In order to determine 03 is is therefore advantageous 
to consider a number of modes. In addition to the dif- 
ferent final states such as those considered above, we can 
also replace K — > K* and D — > D*. Of course if we 
do both at once so that both sides have J ^ then we 
need to do-a more complicated angular analysis as con- 
sidered byl!3. Note that here 03 is common but each case 
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Figure 2: The curves in the </>3 — 6 plane using the parameters 
considered in 10 are shown for the modes K + ir~ (solid curve); 
Kstt (short dashed curve); K + p~ (long dashed curve); K + a^ 
(dash-dot curve); K a p° (dash-dot-dot curve) and K* + x~ (dash- 
dash-dot curve). The inner edge of the shaded region corresponds 
to the 68% CL for TV = 10 8 while the outer edge corresponds to 
the 95% CL. 



has a separate b — axis. To Drive up additional modes 
we can also consider analogous decays where we replace 
the spectator with a d (i.e. B° — > D°K°) or a ud (i.e. 
Aft — * D°A). The point here is that we may be justified 
in putting these cases on a common b axis. 

Because the main point of combining multiple modes 
is to overcome the lack of knowledge of b, B° — > K°D 
where the D is decays to a state such as K + ir~ may be a 
particularly important mode to use in this way since the 
dominant contribution is proportional to b while a is-jcolor 
suppressed in this channel 113. As emphisized byE-3, the 
complimentary case where B~ — > K~ D° and D° decays 
to a CP eigenstate (e.g. n + ir~) and so the a channel is 
much larger than the b channel also gives the same kind 
of information since the amount of interference evident 
in the system determines b without a strong dependence 
on 3 . 

An additional source of constraints that can be help- 
ful may be obtained from a charm factory data which 
can constrain the strong phase differences between D° 
and D° decays as well as S.' V P definitive information con- 
cerning DD oscillations E£EZI. 

Note that a number of the modes we consider are 
instances of three body final states. For a single three 
body mode (e.g. K + ir~ir°) we can consider each of the 
points on the dalitz plot as a having a separate strong 
phase so clearly in principle there is enough information 
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Figure 3: The locus of points on a dalitz plot for the final state 
K + n~ it where Q m i n = Q for 03 = 60° and an overall strong 
phase difference of 0° (solid line) and C = 60° (dashed line). Here 
s = (p n - +Pk+) 2 and t = + p^o) 2 ■ 



to determine_03 . One can thus fit the data to a resonance 
model as inE3 together with the overall strong and weak 
phase differences. 

In this talk I would like to emphasize a different 
method of analysis based on the saturation of eq. ([!]). 
Regarding each point of the Dalitz plot as a separate 
mode, one may find the value of Q m in in Eqn. (Q) for 
each point of the Dalitz plot. Just as in the case of a 
number of discrete modes, the true value of Q must ex- 
ceed all lower bounds. In this case, however, because the 
strong phases due to resonances vary across the Dalitz 
plot, it is likely that the greatest value of Q m in is in fact 
equal to Q. 

In Fig. [| I show a map "magic" points where Q = 
Qmin on the Dalitz -plot for the case of D a — > AT + 7r~-n£ 
using the model ofE£l that uses the data from E687E3 
as input together with SU(3) to give the DCS channel. 
Here I have taken 03 = 60° with an overall strong phase 
difference of 0° for the solid curve and 60° for the dashed 
curve. 

4 The Influence of DD mixing 

In the discussion so far we have explicitly assumed that 
D°D° was negligible. In particular, since we often take 
advantage of interferences involving DCS decays which 
are 0(1%) of the interfering CA decay, the total proba- 
bility of mixing must be less than 0(1%) for this case to 
remain valid. In fact, it has been suggested that the stan- 
dard model may cause DD mixing at about this levelllfl. 
In fact, it can be shown E3 that the changes to d and d 
from such mixing will be O(10%) which leads to an inher- 
ent error in the determineation of 03 of ~ 10° — 15° (p). 
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In order to overcome this possible systematic error, 
there are two approaches: 

1. Using information on the the time between the B~ 
decay and the subsequent D° decay, then the effects 
of possible mixing can be eliminated. 

2. If the parameters of DD mixing are known inde- 
pendently, then they can be taken into account in 
interpreting the time integrated data 

Indeed, if the mixing parameters and time dependent 
data is available, then one can in principle extract ^3 
from just one mode though most likely, the time depen- 
dence in the decay is too weak to make this a useful 
method. 

Here I will emphasise the fact that if we have time de- 
pendent data, it is particularly simple to separate out the 
contributions of mixing and thus proceed with the anal- 
ysis as in the absence of mixing at some cost in statistics. 

The key point is that for D° the decay time is much 
shorter than the oscillation time and therefore it is valid 
to write 

^d(X)^(d (X) + d 1 (X)r)e- T 
dr 

^-d(X) » (d (X) + Ji(X)r)e- T (3) 
dr 

where r = tTi>. Thus do(X) and do(X) would be the 
branching ratios absent mixing so if we extract them from 
the time dependence we may proceed as if there were no 
mixing. 

This can be accomplished through weighting the data 
with wq (t) = 2 — t so that 

POO POO 

d = / [d(r)tflo(r)]dr; d = / [d(r)^o(r)]dr (4) 
Jo Jo 

Using this method more data would be required to 
obtain the same statistical results as in the unmixed 
case. In the unmixed case where d\ — one could ob- 
tain do more effectively by taking the time integrated 
rate. Thus in the unmixed case if a measurement of 
do is based on n events, the uncertainty in do is given 
by: (Ado) 2 /do = 1/n. In the mixed case, using eqn. (||) 
the uncertainty is (Ad Q ) 2 /dl — (2/n) (1 + di/d ) 2 . Since 
di ~ O(do/10), this means that roughly twice the data 
is needed to have the same statistical power as in the 
unmixed case. In order to gauge the precision of time 
measurement required, we can smear out the distribu- 
tion in eq. (W) with a Gaussian resolution function of the 

(T-r'f 

form r(r, r ) oc e ^ where r is the actual time of 
the decay, t' is the measured time of the decay and a is 



the resolution (all in units of 1/T]j). Since r is symmet- 
ric under r r', the fact that w is linear in r implies 
eq. (Q) will still be true for r' but now the error is: 

(Ad ) 2 /d 2 = (2/n)(l + a 2 ) (1 + dx/do) 2 (5) 

As can be seen, the number of events required is not 
adversely effected if a < 1/Td but will be significantly 
degraded otherwise. 

5 Conclusions 

In conclusion, in the system direct CP violation in 
the decay D°K :iz may provide a means of determining 
</>3 with 10 s — 1 9 B mesons. The key is to observe the 
correct D° or combination of D° decay modes. If large 
CP violation is seen in any one mode, it may establish a 
lower bound on sin 2 </>3 while data from multiple modes 
or three body modes can be used to determine sin 2 ^3. 
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